Testing the theory of QGP-induced energy loss at RHIC and the LHC 
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We compare an analytic model of jet quenching, based on the GLV non-Abelian energy loss 
formalism, to numerical results for the centrality dependent suppression of hadron cross sections in 
Au+Au and Cu+Cu collisions at RHIC. Simulations of neutral pion quenching versus the size of 
the colliding nuclear system are presented to high transverse momentum pt- At low and moderate 
Pt, we study the contribution of medium-induced gluon bremsstrahlung to single inclusive hadron 
production. In Pb+Pb collisions at the LHC, the redistribution of the lost energy is shown to play 
a critical role in yielding nuclear suppression that does not violate the participant scaling limit. 
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I. INTRODUCTION 

Strong suppression of single inclusive pions and 
charged hadrons at large transverse momentum, as large 
as Pt = 20 GeV 0, is arguably one of the most 
fascinating phenomena from the heavy ion program at 
the Relativistic Heavy Ion Collider (RHIC). It signifies 
the transition from soft, collective and strongly model- 
dependent physics to the high pT or Et production of 
particles and jets that is well understood in terms of the 
perturbative Quantum Chromodynamics (pQCD) factor- 
ization approach. Calibrated hard probes can thus be 
used to sample the properties of the medium created in 
collisions of heavy nuclei 0. With capabilities to de- 
tect Pt ~ 50 GeV pions and Et ~ 500 GeV jets with 
good statistics, experiments at the Large Hadron Col- 
lider (LHC) will be able to critically test such perturba- 
tive calculations of hadron and jet modification in the 
quark-gluon plasma (QGP) at a new energy frontier. 

Following the discovery of jet quenching in Au-|-Au 
collisions and its verification through d+Au measure- 
ments at RHIC, heavy ion theory has emphasized the 
need for a systematic study of the energy and system size 
dependence of leading particle attenuation 0, Q . Previ- 
ous measurements have been limited to transverse mo- 
menta Pt < 10 GeV and Au-fAu collisions at RHIC 0- 
More recent results up to PT - 20 GeV d Q in both 
Au-|-Au and Cu-|-Cu reactions were shown to be com- 
patible with several theoretical estimates, emphasizing 
the final state QGP-induced suppression of jets [EIEIM,- 
In this Letter, we compare for the first time analytic Q 
and numerical models of jet absorption in order to 
establish the role of sub-leading effects, such as the run- 
ning of the strong coupling constant as with the Debye 
screening scale /i ~ gT in the plasma. We present our 
predictions for tt" quenching versus centrality, pT and 
^/s, based on the Gyulassy-Levai- Vitev (GLV) approach 
to the medium-induced non-Abelian energy loss 10'| . We 
note that redistribution of the lost energy in small and 
moderate pT hadrons plays a significant role in the tran- 
sition from high pT suppression to low pT enhancement 
of the back-to back two particle correlations jllj . In this 



Letter we give results for single inclusive pions and iden- 
tify the range of transverse momenta at RHIC and the 
LHC where the gluon feedback is important. 

In Section II we present a simple analytic model for 
the system size, iVparti dependence of jet quenching in 
nucleus-nucleus collisions. Section HI contains select in- 
termediate results from the GLV theory, including the 
evaluation of the mean gluon number (Ng) and frac- 
tional energy loss (AE/E). The derivation of the radia- 
tive gluon contribution to small and moderate pT hadron 
production is given in Section IV. Section V presents the 
calculated suppression of 7r° production in Cu-|-Cu and 
Au-|-Au collisions at RHIC and Pb-|-Pb collisions at the 
LHC. Conclusions are given in Section VI. 



II. ANALYTIC MODEL OF JET QUENCHING 

While tomographic determination of the properties of 
the medium created in nucleus- nucleus collisions can only 
be achieved through detailed numerical simulations, it is 
useful to define a simplified analytic model which incor- 
porates the essential features of jet quenching calcula- 
tions. We first review the approach formulated in |^ and 
discuss its advantages and limitations. 

The final state medium-induced energy loss occurs af- 
ter the hard partonic scattering ab — *■ cd and before the 
fragmentation of the parent parton into a jet of colorless 
hadrons. We consider for simplicity high pT > 5 GeV 
hadron production at RHIC and will show later that the 
contribution of the induced bremsstrahlung in this region 
is small. There are two equivalent ways of implementing 
energy loss in the perturbative QCD hadron production 
formalism. The first one associates AE with the kine- 
matic modification of the momentum fraction z = p'^ /p+ 
in the fragmentation function Di^/^(z), leaving the hard 
parton production cross section da'^ / dydPpT^ unmodified. 
The second approach reduces the jet cross section in the 
presence of the medium but leaves Di^/^(z) unaltered. It 
can conveniently be implemented in an analytic model 
of QGP-induced leading hadron suppression. We take 
the underlying parton production cross section to be of 
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power law type 
da'' 



A 



, if PT, >po , (1) 



where n = n {y,pT^, y/s) and po ~ 2.5 GeV at RHIC. 
In a finite px range, fixed rapidity y and center of mass 
energy ^/s, a constant n = {n {y,pT^, ^/s)) is a good 
approximation. At RHIC, in the region of 5 GeV < 
PTe < 10 GeV, the spectra scale roughly as Uq — 7, 
Ug = 8.4. Fragmentation functions are convoluted with 
the partonic cross section as follows 

da'' 
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In Eq. ^ Zmin = Pt/pTc max and the subsequent ap- 
proximation is most reliable when Zmin <C 1- It should 
be noted that (z) will also depend on the partonic and 
hadronic species. 

The second input to the analytic model comes from 
the radiative energy loss formalism Hj]. To understand 
the system size dependence of jet quenching, we use the 
approximate GLV formula which relates AE to the size 
and the soft parton rapidity density of the medium. For 
(1-|-1)D Bjorken expansion, in the limit of large parton 
energy 2E/n'^L > 1, we find O 
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In Eq. -/j is the jet path length in the medium and A±^ 
is the transverse area. Cr — 4/3 (3) for quarks (gluons), 
respectively, is the quadratic Casimir in the fundamental 
(adjoint) representation of SU(3). Numerical simulations 
of AE/E clearly indicate a weaker dependence of the 
fractional energy loss on the jet energy than given in 
Eq. ©. 

The key to understanding the dependence of jet 
quenching on the heavy ion species and centrality is the 
effective atomic mass number, Aes, or the number of par- 
ticipants, A'part, dependence of the characteristic plasma 
parameters in Eq. © j3|, 



dN^/dy cx dN^/dy cx A^s cx iVpart , 
A^ cx A%^ 
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Therefore, the fractional energy loss scales approximately 
as 



e = AE/E cx A 



2/3 



2/3 



part 



(6) 



up to logarithmic corrections from Eq. (PJ- If a parton 
loses this momentum fraction e during its propagation in 



the medium to escape with momentum p^"'^"'^'*, immedi- 
ately after the hard collision pr^ = p^"'^"'^'*/ (1 — e). Not- 
ing the additional Jacobian \d'^p'!^'^^'''' / d^pTj\ = (1 — e)^, 
we find for the quenched hadronic spectrum per elemen- 
tary NN collision 
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In Eq. O Ecff is the average over all parton species and 
accounts for the color charge, geometry and multi-gluon 
fluctuations. From this result we can easily derive the 
system size dependence of the nuclear modification factor 



Raa = 



a;i dN'iJdyd^pT 
Naa col da^/dyd^pT 
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In Eq. © k/{n — 2) is the proportionality coefficient in 
Eq. which depends on the microscopic properties of 
the medium but not on its size. With n 3> 1 and ex- 
perimentally measured and theoretically calculated sup- 
pression, ~ 5 fold in central Au-|-Au collisions at RHIC, 

Ecff = -^part^/('^ ~ 2) is Small. We thus predict that 
the logarithm of nuclear suppression, Eq. (jSJ, has simple 
power law dependence on the system size 



In Raa = 



-IcN"^^^ 

'>'^*part 
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where the leading correction goes as k'^N^l^^/2{n — 2). 

We focus on the most central nuclear collisions and 
use impact parameters b = 1 — 3 fm depending on the 
atomic mass A. An optical Glauber model calculation is 

2/3 

used to evaluate A^part with cr™ = 42 mb and a Wood- 
Saxon nuclear density with results given in Table The 
analytic prediction for the system size dependence of jet 
quenching is shown in Fig. ^ It is fixed by the mag- 
nitude of the suppression established in central Au-|-Au 
collisions and consistent with existing simulations ^3 ■ 
From this analysis we expect a factor ~ 2 suppression in 
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TABLE I: Summary of the relevant quenching parameter 
-^part (rounded) at a fixed impact parameter b for select heavy 
ion species. 
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FIG. 1: Top panel: the predicted linear dependence of 
jet quenching in natural variables, In Raa versus N^i'i , for 
central collisions of '"'Be, 
and ^^*U. Bottom panel 
to PHENIX and STAR experimental data in Au+Au and 
Cu+Cu collisions. 
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central Cu+Cu collisions. Comparison to the preliminary 
PHENIX pT > 7 GeV data [J and STAR pr > 6 GeV 
data 2] is also shown in the bottom panel of Fig. with 
good agreement within the experimental uncertainties. 
Dashed and solid lines illustrate the difference between 
Eqs. (IHI and @ when normalized to the same suppres- 
sion in central Au+Au collisions. 

The main advantage of the GLV analytic model is the 
ability to provide guidance on the magnitude of the QGP- 
induced jet quenching versus centrality and address a 
large body of experimental data. Its limitations include 
a fixed coupling constant as , implementation of only the 
mean, though suitably reduced, energy loss AE and the 
inability to incorporate additional nuclear effects, such 
as the Cronin multiple scattering j5(] and nuclear shadow- 
ing . It also relies on a reference numerical calculation 
in central A+A collisions in the same pr and y range as 
well as ^/s The deviation of dN^ /dy from the exact 



participant scaling in Eq. Q may lead to less quench- 
ing and improved agreement with the data in peripheral 
reactions, but is here neglected. 



III. NUMERICAL EVALUATION OF THE 
QGP-INDUCED ENERGY LOSS 

The solution for the differential in energy (uj) and 
transverse momentum (k) spectrum of medium induced 
gluon radiation has been obtained order by order in the 
correlations between the multiple scattering centers in 
nuclear matter using the reaction operator approach [lol | 
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Here are the momentum transfers from the medium, 
distributed according to a normalized elastic differential 
cross section aei{i)^^d(Jei(i) / d'^cii, and Az^ — Z}~ — Zk-i 
are the separations of the subsequent scattering centers. 
In Eq. (|10() the color current propagators and inverse 
formation times are denoted by 
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Numerical results have been obtained to third order in 
the opacity, (L/Xg), with all correlations of up to four 
scattering centers, including the initial hard interaction. 
To speed up the evaluation of the squared amplitudes, 
the oscillating quantum coherence phases between the 
points of interaction were converted to Lorentzians using 
an exponentially falling geometry with no sharp edges. 
The requirement that the mean locations Zk, k — \ ■ ■ - n 
of n scatterers are the same as the one for a uniform 
soft parton distribution fixes the parameter L,, of such 
geometry 10]: 



{zo - Zk) 



L,. 



(12) 



Bjorken (1+1)D expansion of the plasma is accounted for 
as follows: 



p{zk) = p{zo) 



£0 



/ ^ 1/3 /X 1/3 

fi{zk) = ^J-{zo) i^j , Ag(zfc) = Ag(zo) f ^ j ,(13) 
and the kinematic constraints, 



^i/Q <x = k+/E+ w uj/E < 1 
M 5: |k| < \/Q^ mm{x, 1 — a;) , 



(14) 



have been incorporated for consistency with our previous 
work |l2l|. It was recently shown that for physical on- 
shell final state gluons the medium induced radiation is 
infrared and coUinear safe . This allows relaxation of 
the fi/Q < X, fi < |k| constraints in the future, though 
it should be noted that the Debye screening scale still 
controls the small k and uj cancellation between the single 
and double Born diagrams in the opacity expansion [lOj | . 

To evaluate the effect of multiple gluon emission and 
arrive at a probabilistic distribution P{e) for the frac- 
tional energy loss e = AE/E = X^iLi ^' ^i/E, we 
are motivated by an independent Poisson gluon emission 
ansatz, but incorporate kinematic constraints [T^ : 



Pnie) 



I r de' P^^,{e-e')^ie' ^ij/E) . (15) 
n Jq de' 



We normalize this probability density to unity and 
Eq. H15|l ensures that at every step energy is conserved. 
As a consequence, for small jet energies and large AE/E 
the gluon distribution is distinctly non-Poisson. We eval- 
uate the mean energy loss as follows: 



(16) 



In Ref. .12j we considered jet production following the 
binary collision density Ta^(6) in central Au-|-Au re- 
actions at ^/s — 200 GeV. In an elementary hard in- 
teraction inclusive jets are distributed uniformly in az- 
imuth relative to the reaction plane. We calculated 
((A£'))goom. using the hue integral, Eq. (^UJ, through 
the (1-|-1)D Bjorken expanding medium density by cor- 
rectly weighing the amount of lost energy with the jet 
production rate. Cylindrical geometry with radius L = 
6 fm that gives the same mean energy loss ( AE ) for 
uniform initial soft parton rapidity density was then con- 
strained. Using Eqs. I@J and ^ we can determine the 
effective length L, transverse area A± = ttL^, gluon ra- 
pidity density dN^/dy and effective atomic mass A^s 
for shadowing applications 14] in interacting heavy ion 
systems of different size. Table ^ summarizes the pa- 
rameters used in our calculation of central, mid-central 
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FIG. 2: Left panel: mean energy loss for Au+Au and Cu+Cu 
soft parton densities at RHIC corresponding to centralities 
given in Table ^ and Pb+Pb soft parton densities at the 
LHC corresponding to central collisions versus the jet energy. 
Right panel: mean gluon number versus Ej^t for the same 
heavy ion systems. 



and peripheral Au+Au and Cu+Cu collisions at ^/s = 
200 GeV at RHIC. For central collisions with ^/s — 
5.5 TeV at the LHC, we use dN^/dy = 2000, 3000 and 
4000 to test the sensitivity of jet quenching to the QGP 
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TABLE II: Estimated dN^/dy, L and AeS versus A'^part for 
central, semi-central and peripheral Au+Au (top table) and 
Cu+Cu (bottom table) at RHIC. 



properties. 

The calculated fractional energy loss and mean gluon 
number for quark and gluon jets versus their energy for 
the centralities and densities discussed above are shown 
in the left and right panels of Fig. |21 Only in the 
limit Ejet oo, AE/E does the energy loss for 
quarks and gluons approach the naive ratio AE^ /AE'^ — 
Ca/Cf = 9/4. For large fractional energy losses this ra- 
tio is determined by the AE < E constraint. In this 
regime no simple scaling arguments related to the prop- 
erties of the dense nuclear matter and the color charge 
are applicable. It should be noted that at high Ej^t the 
fractional energy loss is not large even at the LHC. 

In our calculation the strong coupling constant as is 
not used as a free parameter but evaluated at the typical 
scale fi in the elastic scattering cross section (Jei{i) and 
the gluon mean free path Xg{i) — 1 / a ei{i) p{i) ■ At the ra- 
diation vertex, as(k^) is also sensitive to the increase of 
the temperature or density of the medium, p (xT^, with 
the increase of dN^ /dy at a fixed transverse area A±. 
The effects described here lead to sub-linear dependence 
of the energy loss on dN^ /dy. For the LHC example, 
given in Fig. |21 this can be a 50% correction for large 
Ejet and even more significant at low Ejet when com- 
pared to Eq. |j2Jl. Conversely, at low parton densities in 
peripheral collisions the energy loss will be larger than 
naively expected. 

Another important point, seen in the right panel of 
Fig. 121 is that except for very low jet energies the mean 
gluon number (Ng) is not small (<C 1) and the probability 
of not radiating gluons, Pq = exp(— (iVJ), is never large. 
Full results for P{e) were shown in jMIQ- ^^'^ that 
the probability density does not approximate aS{e) -I- (1 — 
a)S{l — e), a < 1. The latter ansatz yields RAAiPr) = a 
independent on the collision energy or px- Instead, P{e) 
is much more uniformly distributed in the interval [0, 1] 
and our calculations retain sensitivity to the local slope 
and the parton species contribution to the differential 
inclusive hadron production cross section, see Eq. (jSJ. 

We finally note that multi-gluon fluctuations, given by 
Eq. (|15|l . reduce the jet quenching effect relative to the 
application of the mean ( AE/E ) shown in Fig. 1^ [Tsj. 
This can be seen by comparing the Ccff ~ 0.2 in central 
Au-|-Au reactions, obtained from Eq. ||SJ), to the frac- 
tional energy loss for 10 GeV quark jets, Ccff < ( AE/E ). 
To investigate the scaling of energy loss with iVpart we 
select quark jets of Ej^t — 10, 20 and 30 GeV at midra- 
pidity at RHIC and A^part = 9,21,55,103,167 and 328. 
These cover fractional energy losses < ( AE/E ) < 0.4 
with numerical results shown in Fig. |3| Power law fits, 
also plotted, give exponents n = 0.60 — 0.63 that arc not 
very different from the naive n = 2/3 expectation from 
Eq. lO that was used in Fig. ^ We conclude that the 
deviation between the calculated energy loss with kine- 
matic constraints and running strong coupling constants 
and its asymptotic flxed as analytic behavior is smaller 
when the variation of the energy loss is associated with 
a change in the system size rather than a large change in 
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FIG. 3; A'^part dependence of the fractional energy loss AE/E 
for E-jet = 10, 20 and 30 GeV quark jets in = 200 GeV 
nucleus-nucleus collisions. Power law fits yield exponents n — 
0.61,0.60 and 0.63, respectively. 



the density of the medium alone. 



IV. NUCLEAR EFFECTS ON INCLUSIVE 
HADRON PRODUCTION 

The lowest order perturbative QCD cross section for 
single inclusive hadrons in nucleon-nucleon (NN) reac- 
tions, including non- vanishing transverse momentum ka_h 
distributions of the incoming partons, is given by 



dyd^PT 



X /(kQ)/(kf,) (f>a/N{Xa,Pf)(t>b/N{xb, M/) 

1 da"''-""^ 
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TTZr 



dt 



-Dh/c{Zc,f^f) 



(17) 



Here Zc = Pt/pt,, Xa,b = Pa,b/ 
notation <t>a.b/N{xb, P'f) are 



In 



the parton distribution 
functions (PDFs), Di^/f.{zc, p.f) are the fragmentation 
functions (EEs) and we have chosen the factoriza- 
tion/fragmentation, and renormalization scales /z/ = 
Mr = PTc- The distribution of non-zero transverse mo- 
menta of the incoming partons is parametrized as follows 
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^{klb. 



■eM-Kb/iklb))- 



(18) 



The physical requirement for hard partonic scattering is 
ensured by ka.b < Xa^b\fs and K = 1.5 is a phenomeno- 
logical K-factor at RHIC. For further details see P, 5. 
Comparison to the PHENIX measurement of 7r° produc- 
tion in -/i = 200 GeV p-l-p collisions at RHIC is shown 
in the insert of Fig. ^ 

Nuclear effects can be incorporated in the pQCD for- 
malism, Eq. (|17|) . and fall in two categories: medium- 
induced kinematic modifications to the perturbative for- 
mulas and possibly universal modifications to the PDFs 
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FIG. 4: The predicted invariant multiplicity distribution of 
neutral pions in central Au+Au collisions at ^/s = 200 GeV 
for medium density dN^/dy = 800 — 1175 and TauAu = 
23 m6"\ The same calculation for Cu+Cu collisions for 
medium density dN^/dy = 255-370 and TcuCu = 4.5 mfe"^ 
The insert shows the cross section for 7r° production in p+p 
collisions to LO pQCD. Data is from PHENIX 0. 



and FFs. An example of the latter are parameterizations 
of nuclear shadowing Sq,g/A{x, fif) Il4j included in this 
calculation via 



A 



A 



N 



+ -JSq,a/Aix,^if)(|}q,g/nix,^if)] .(19) 



In future work we will combine dynamical calculations of 
coherent nuclear enhanced power corrections with other 
elastic and inelastic effects in nuclear matter [l^ . Trans- 
verse momentum broadening of the incoming partons, 
leads to enhancement of the differential particle distribu- 
tions at Pt ~ few GeV 5] and can be accounted for in 
Eqs. (HZj), (0 as follows: 



{kit) 



{ka,b) vac ' 



{k-a.hlmea 



(20) 



Here {kl^^)jned = i'^fJ-'^La^b/ K.b)^ and the typical momen- 
tum transfers squared /i^, mean free paths Xa^b and par- 
ton propagation lengths La^b refer to cold nuclear mat- 
ter gig. 

In the QGP, final state energy loss is the dominant 
effect that alters the single and double inclusive hadron 
production cross sections |^ 0| . Application to the at- 
tenuation of leading hadrons as a kinematic modification 
of the momentum fraction z in the FFs Dii/c{z) is con- 
sidered standard H, ■ The redistribution of the lost 
energy in soft and moderate pt hadrons was only recently 
derived in the pQCD approach, first for back-to-back 
two particle correlations [111. The established dramatic 
transition from the high px factor of four suppression 



{Raa ~ 0.25) to the low pT factor of two enhancement 
{Raa ^ 2) makes it imperative to study this effect for sin- 
gle inclusive particle production. One possibility is that 
the fraction of the hadrons from the bremsstrahlung glu- 
ons is negligible or small over the full accessible px range. 
At the other extreme, a very large fraction may compro- 
mise the current jet quenching phenomenology, leading 
to Raa ^ 1 at moderate transverse momenta even in 
dense matter. 

With this motivation, we first give results for the mod- 
ification of inclusive hadron production from final state 
radiative energy loss. It can be represented as: 

Here P(e) is calculated from Eqs. IpH, and dNg/de 
is the distribution of the average gluons versus e — uo/E 
so that 



de -j-ie) ^ (Ng) , I de e^{e) 



de 



de 



(22) 

It is easy to verify the momentum sum rule for all 
hadronic fragments from the attenuated jet and the ra- 
diative gluons. With appropriate changes of variables 

^ dz zDh/c{z) =^ 

f\ dN9 . ^ r dz Z ^ (Z\ 



= 1 - e + e = 1 



(23) 



Figure shows the invariant tt'^ multiplicity in 
central Au+Au and Cu+Cu reactions. Data is from 
PHENIX 0]. At high pr, comparisons been the jet 
quenching theory and the experimental measurement can 
(and should) also be made for the differential cross sec- 
tions. At low pT^ a deviation is present due to the fixed 
order baseline pQCD calculation and QGP effects are 
more accurately studied via Raa{pt)- 



QUENCHING OF INCLUSIVE PIONS AT 
RHIC AND THE LHC 



Having evaluated the energy loss of quark and gluon 
jets in the QGP media specified in Table ^ we calcu- 
late the quenched pion spectra using Eqs. lfT7|l . lfTO|l . 
(Pn|l and (|?T|l . Figure [S] shows Raa{pt) for central, 
semi-central and peripheral collisions. The predictions 
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FIG. 5: Top panel: nuclear modification factor Raa for 
Au+Au collisions at -^s = 200 GeV versus pr and central- 
ity. Bottom panel: similar calculation for Cu+Cu at _v/s = 
200 GeV at RHIC. Preliminary data is from PHENIX fO, and 
analytic estimates from Fig. Q are also shown. 
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FIG. 6: Suppression of tt" production in central Au+Au 
collisions at the LHC as a function of the parton rapidity 
density. Insert shows the baseline p+p 7r° cross section at 
= 200 GeV and = 5.5 TeV. 



in Cu+Cu reactions are for a constant suppression fac- 
tor, as in Au+Au, at higli pT- Preliminary PHENIX 
data, first compared to this theory in Ref. are also 



included, ft should be noted that even in peripheral reac- 
tions there can be noticeable particle attenuation. This 
is larger than the analytic estimates due to sub-leading 
effects of the medium density on the parton energy loss, 
discussed in the previous Section, see Fig. Whether 
such effects are observable or compensated by the non- 
uniform QGP density in the transverse plane is subject 
to experimental verification. Finally, we make the im- 
portant observation that for similar densities and system 
sizes, for example dN^ /dy = 410 in mid-central Au+Au 
and dN^ /dy = 370 in central Cu+Cu, the magnitude of 
the predicted pion suppression is similar. 

Pb+Pb collisions at the LHC represent the future en- 
ergy frontier of QGP studies in heavy ion reactions. We 
have explored the sensitivity of Raa{pt) to the par- 
ton rapidity density in central nuclear reactions with 
dNS/dy ~ 2000,3000 and 4000. In J a seven-fold in- 
crease of the medium density in going form RHIC to the 
LHC was assumed. In this work we adhere to a more 
modest two- to four-fold increase of the soft hadron rapid- 
ity density and emphasize that future measurements of 
jet quenching must be correlated to dN^ /dy ~ dN^ /dy ~ 
{3/2)dN''^/dy 0, H [13 to verify the consistency of the 
phenomenological results. 

From Figs. [3 and [7| we conclude that at low and 
moderate px, jet suppression at the LHC is larger than 
at RHIC. However, at high pT > 30 — 50 GeV this order- 
ing is reversed. The physics reason for our result is that 
the fractional energy loss (AE/E) of fOO GeV quark 
and gluon jets at the LHC is not very different from that 
of 25 GeV jets at RHIC, see Fig. El In addition, the 
power law dependence of the hadronic and the underly- 
ing partonic spectrum n (2/,PTe7 a/s) changes (decreases) 
from = 200 GeV to = 5.5 TeV. This is shown in 
the insert of Fig. and affects the magnitude of calcu- 
lated nuclear suppression, Eq. LHC will soon provide 
an extended px range to test the correlation of RAAiPr) 
with the stiffness of the differential particle spectra. 

To assess the importance of the gluon feedback term in 
Eq. (|21|l . we extend, in Fig.[7| the calculation of Raa{pt) 
for central Au+Au and Pb+Pb collisions to low and mod- 
erate transverse momenta. At RHIC the redistribution of 
the lost energy leads to small, ~ 25%, modification of the 
neutral pion cross section in the pt ~ 1 — 2 GeV range 
and modest improvement in the theoretical description 
of that data. In contrast, at the LHC the fragmenta- 
tion of medium-induced gluons is a much more signifi- 
cant > 100% correction to the low and moderate pr Tr*^ 
production rate. 

The need for the more consistent treatment of jet 
energy loss, Eq. H21|l . is also illustrated by compar- 
ing RAAipr) in Fig. [7| to the participant scaling ratio: 
(-^part/2)/-/Vcoi 0|. At high Pt there is no lower limit 
on the quenching of jets. At low pT the total available en- 
ergy of the collision, ~ {Npa,rt/2)y/s, suggests participant 
scaling of bulk particle production confirmed by hydro- 
dynamic calculations Previous estimates of leading 
particle suppression at low and moderate pr have vio- 
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FIG. 7: Top panel: nuclear modification factor Raa in 0-10% 
central Au+Au collisions at moderate pr with (solid line) and 
without (dashed line) gluon feedback, dN^ /dy = 1175. Pre- 
liminary data is from PHENIX Bottom panel: nuclear 
suppression at moderate pT at the LHC y/s = 5.5 TeV. Cen- 
tral Pb-|-Pb collisions with (solid line) and without (dashed 
line) gluon feedback are shown, dN^/dy ~ 2000, 3000, 4000. 



lated this limit [a,l3- Fig. [7| shows that the gluon feed- 
back can ensure numerically RAAiPx) > (-^part/2)/-/Vcoi 
at mid-rapidity at the LHC for the densities consid- 
ered here and is important everywhere in the region 
PT < 15 GeV. 



centrality and transverse momentum. Pb-|-Pb reactions 
at — 5.5 TeV at the LHC were also discussed and 
our calculations at mid-rapidity accounted for Cronin 
multiple scattering 0, nuclear shadowing ^Jj and final 
state radiative energy loss in the quark-gluon plasma 0| . 
Elastic energy loss was not considered here, since its ef- 
fects are still under debate ^3 relative to the attenua- 
tion of jets via gluon bremsstrahlung. We compared our 
numerical results to a simplified analytic model for cen- 
trality dependent jet quenching and showed that there 
are non-negligible corrections in the evaluation and im- 
plementation of radiative energy loss related to the tem- 
perature or /i dependence of the strong coupling constant. 
While the dependence of the observable hadron suppres- 
sion on dN^/dy was shown to be sub-linear, this calcu- 
lation retains sensitivity to both the properties of the 
medium and the underlying perturbative baseline cross 
sections. 

At low and moderate transverse momenta we derived 
the contribution to single inclusive pion production from 
the fragmentation of medium-induced gluons. At RHIC 
we found this effect to be a modest, < 25%, correction. 
Our result should be contrasted with the case of back-to- 
back di-hadron correlations where the redistribution of 
the lost energy controls the QGP-induced transition from 
suppression to enhancement of large angle inclusive two 
particle production ,ll|. At the LHC, however, even in 
inclusive one pion calculations, gluon feedback is shown 
to alter the < 15 GeV cross section in central Pb-|-Pb 
reactions by as much as a factor of two. In summary, 
results reported in this Letter not only provide a more 
consistent theoretical framework to treat the effects of 
medium- induced gluon bremsstrahlung but also rectify 
the over- quenching of jets at low px in the limit of large 
fractional energy loss 



jets 1 

m 



VI. CONCLUSIONS 

In this Letter we presented predictions for the nu- 
clear modification of inclusive neutral pion production in 
Au-|-Au and Cu+Cu collisions at ^/s = 200 GeV versus 
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